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ABSTRACT:  Analysis  of  ab  initio  wave  functions  shows  that  the  interaction  between 
halogen  atoms  (F,  Cl,  and  Br)  and  Ag  clusters  is  ionic,  and  that  the  halogen  ionicity  is 
essentially  —1.  The  interaction  and  the  bonding  arise,  almost  entirely,  from  the 
Coulomb  attraction  between  the  charged  halogen  and  the  metal  and  from  polarization 
of  the  two  subunits.  Large  shifts  in  equilibrium  bond  distances  and  vibrational 
energies  are  caused  by  an  external  electric  field.  These  changes  arise  form  a  dominant 
Stark  effect. 

‘Present  Address:  Dipartimento  di  Chimica  Inorganica  e  Metallorganica,  Universita  di 
Milano,  20133  Milano,  Italy 
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The  nature  of  the  bonding  between  halogen  atoms  and  silver  is  a  subject  of 
considerable  interest  in  surface  science  and  electrochemistry  [1-4],  and  the  nature  of 
the  bonding  -interaction,  if  more  ionic  or  more  covalent,  is  still  a  matter  of  controversy. 

We  studied  the  interaction  between  F,  Cl,  and  Br  atoms  and  an  Ag4  cluster 
modeling  the  three-hollow  site  of  the  Ag(lll)  surface  using  ab  initio  Hartree-Fock 
wavefunctions.  AtTPntipn  has  been  fccueed  on  the  character  of  the  bonding  and  on  the 
importance  of  the  polarization  efieas  for  the  dipole  moment  changes  resulting  from  the 
bond  formation.  In  addition  we  have  considered  the  effect  of  an  external  electric  field 
on  the  metal-adsorbate  properties. 

A  major  concern  in  adsorbate  cluster-model  studies  is  the  dependence  of  the 
properties  on  the  number  of  substrate  atoms  in  the  cluster.  However,  the  analysis  of 
the  interaction  of  F  atoms  with  Agg,  Ag7,  and  Agjg  clusters  confirms  the  basic  features 
of  the  bonding  derived  for  Ag4  [5].  Here  we  discuss  only  the  results  for  Ag4,  but  we 
believe  that  our  conclusions  about  the  character  of  the  bonding  will  be  valid  for  more 
general  cluster  geometries. 

The  molecular  wave  functions  have  been  determined  using  flexible  basis  sets  of 
contracted  Gaussian  t>pe  orbitals.  In  order  to  facilitate  the  calculations,  the  core 
electrons  of  Ag  have  been  replaced  by  Effective  Core  Potentials  [6]  including  explicitly 
in  the  valence  only  the  4s^  4p^  4d^^  5s^  electrons.  The  valence  basis  set  for  Ag  is 
contracted  to  [4s4p3d].  The  halogen  atoms  have  been  treated  at  all  electron  level  and 
the  contracted  basis  sets  used  are  [5s4pld],  [5s4p2d],  and  [9s7d3d]  for  F,  Cl,  and  Br, 
respectively  ["].  The  basis  sets  optimized  for  the  neutral  atoms  [7]  were  extended  with 
diffuse  s,  p,  and  d  functions  whose  exponents  have  been  determined  in  order  to  obtain 
a  better  approximation  for  the  SCF  electron  affinity  (EA),  and  for  the  dipole 
polarizabilitv'.  xq,  of  the  respective  anions. 
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The  halogen  atom  to  surface  distance  in  Ag„-X  has  been  varied  along  the  C3y 
axis:  the  distance  of  X  from  the  surface  plane  (r),  the  dissociation  energy  (Dp),  the 
vibrational  frequency  (cug)  and  the  dipole  moment  (fi)  have  been  determined  (Table  1). 

For  distances  close  to  the  equihbrium  Ag-halogen  distance,  the  ground  state 
configuration  is  (lajf  {laif  (3ai)^  (le)"*  (2e)^  (^E);  the  halogen  cores  as  well  as  the 
Ag  Is^  to  4d^®  cores  are  not  explicitly  given.  This  state  arises  from  the  combination  of 
the  Ag4  ground  state  (aj)^  (e)^  with  the  seven  valence  ns  and  np  electrons  of  X  in 
configuration  (a^)^  (e)^.  The  main  bonding  effects  are  in  the  outer  e  shells  and 

can  be  either  covalent  or  ionic.  In  this  latter  case  the  bonding  can  be  considered  as 
originating  from  the  interaction  of  the  ionized  Ag4  (aj)^  (e)*  cluster  and  the  closed 
shell  anion  X  (e)'*.  In  the  covalent  view,  a  filling  bonding  (le)'*  and  a  singly 

occupied  antibonding  (2e)^  .MOs  are  formed. 

To  determine  the  character  of  the  Ag-halogen  bonding  as  covalent  or  ionic  we 
apply  twi  measures.  First,  we  make  use  of  a  projection  operator  [8]  for  the  occupied 
orbitals  of  the  free  halogen  anion  X  .  The  expectation  value  of  the  projection 
operator  for  a  given  spatial  orbital  of  X~,  j  i  m-  is  a  measure  of  the  extent  to 
which  1  m  is  occupied  in  the  total  cluster  wavefunction,  (/'(Ag^X).  In  the  case  that 
this  projection  gives  a  value  close  to  2,  it  is  clear  that  of  X~  is  completely 

occupied  in  ip{A^X)  [8].  For  all  X  and  for  a  wide  range  of  distances  about  the 
minimum  of  the  Ag4-X  energy  curve,  the  projection  of  all  0^  j  is  about  2.  Summing 
over  all  the  orbitals  of  X,  we  find  9.98,  17.95,  and  35.92  electrons  for  F,  Cl,  and  Br, 
respectively,  in  Ag4-X;  these  values  are  for  the  halogen  near  the  minimum.  They 
provide  very  strong  evidence  for  an  ionic  interaction  where  the  silver  cluster  is  Ag4  and 
the  halogen  is  X”. 
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A  second  measure  of  the  bonding  nature  is  based  on  the  slope  of  the  dipole 
moment  curve  as  the  distance  of  X  from  the  cluster  is  varied  [8].  For  an  ionic  bond, 
the  absolute  magnitude  of  the  slope  should  be  —1.  In  fact,  the  non-zero  component  of 
the  dipole  moment  of  an  ionic  molecule  represented  by  two  point  charges  -i-  q  and  — q 
will  be  =  — q  r  and  d/i/dr  =  — q.  Therefore,  for  an  ideal  fully  ionic  molecule  where 
q  =  1,  d/i/dr  =  —1  and  the  curve  is  a  straight  line.  Generally,  on  the  other  hand,  a 
small  slope  is  characteristic  of  a  covalent  bond.  Of  course,  in  a  real  system  the 
polarization  of  the  subunits,  in  our  case  the  halogen  and  the  cluster,  can  lead  to  a 
substantial  change  in  the  absolute  magnitude  of  the  dipole  moment.  In  our  case,  this 
polarization  must  be  quite  important  because  we  find  a  slope  greater  than  1,  see 
Table  1;  it  is  impossible  that  the  ionicity  of  the  halogen  is  greater  than  1. 

The  dipole  moment  curve  as  function  of  the  distance  has  been  determined  by  a 
cubic  polynomial  fitting: 


/i(r)  =  Mo  +  M  ,(r  -  r^)  -I-  MjCr  -  +  .Mjlr  -  re)^ . 


For  each  value  of  r  it  is  possible  to  transform  the  Ag4X  cluster  orbitals  into  those 
which  are  most  like  X,  0i(loc,X),  and  the  orthogonal  complement  of  the  SCF  canonical 
space.  This  corresponding  orbital  transformation  [9,10]  provides  the  transformation 
with  maximally  localized  orbitals.  The  contribution  to  the  total  dipole  moment  of 
0i(loc,X)  to  the  total  cluster  dipole  is  denoted  by  /i(X).  The  coefficients  of  the 
expansions  of/i(Ag4X)  are  reported  in  Table  1.  The  dipole  moment  curve,  for  a 
considerable  distance  about  the  equilibrium  Ag-halogen  bond,  is  a  straight  line  (small 
M2  and  M3)  and  the  contribution  from  the  chemisorbed  X-halogen  to  the  slope  of  the 
curve  is  almost  exactly  —  1  thus  supporting  the  view  of  an  essentially  ionic  bond. 
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In  order  to  understand  the  origin  of  the  Ag^-X  interaction,  we  applied  a 
constrained  space-orbital  variations  CSOV  [II].  In  this  way  it  is  possible  to  determine 
the  individual  contributions  to  the  bonding,  i.e.,  electrostatic  effects,  intraunit 
polarization,  charge  transfer,  and  covalent  bonding.  We  start,  CSOV’^  step  0,  with  the 
superposed  Ag4  and  X~  charge  distributions;  this  gives  the  frozen-orbital  (FO) 
interaction.  Given  the  strong  evidence  that  the  bonding  is  ionic,  we  analyze  the  effect 
of  replacing  the  X  ion  with  a  unit  point  charge,  PC  =  —  1,  placed  at  the  position  of 
the  halogen  nucleus.  In  this  case  step  0  gives  the  pure  electrostatic  attraction  with  no 
repulsive  contribution  originating  from  non-bonding  overlap  of  X  and  .AgJ  charge 
distributions  (Pauli  repulsion).  Then  we  allow  the  .AgJ  orbitals  to  vary  because  of  the 
presence  of  the  X  anion  for  both  the  PC  and  the  extended  representation  of  X  (CSOV 
step  1);  this  allows  the  Ag^  to  polarize  in  the  presence  of  X  .  The  next  step  (step  2)  is 
the  Ag4  orbital  variation  in  the  full  space,  indicating  the  amount  of  charge  transfer 
from  Ag4  to  X”.  In  step  3  the  Ag4  orbitals  are  fixed  and  the  X  orbitals  are  varied  in 
their  own  space.  This  permits  the  X~  unit  to  polarize  as  a  result  of  the  presence  of 
Ag4 .  In  CSOV  step  4  this  variation  occurs  in  the  full  space,  so  that  both  X~  to  Ag4 
donation  and  X”  polarization  are  possible.  In  step  5  the  covalent  mixing  between  the 
valence  e  orbitals,  (le)'^  (X~),  and  (2e)^  (Ag4  ),  is  allowed;  all  other  orbitals  are  fixed  as 
determined  in  step  4.  The  result  of  step  5  is  compared  to  the  full  unconstrained  SCF 
result  and  if  they  are  nearly  the  same  this  means  that  all  important  bonding 
contributions  have  been  considered.  Of  course,  steps  2-6  have  no  meaning  when  X  is 
replaced  by  a  point  charge.  For  each  step,  we  consider  the  interaction  energy  with 
respect  to  separated  Ag4  and  X  fragments,  and  the  dipole  moment In  Table  2  the 
results  of  this  analysis  are  reported  for  Ag4-F  at  r^  =  3.2  near  equilibrium  bond 
distance.  Clearly,  the  most  important  contribution  comes  from  the  electrostatic 
interaction.  The  difference  between  the  PC  and  F  values,  1.5  eV,  represents  the  Pauli 
repulsion  contribution.  The  electrostatic  attraction  is  followed  by  the  cluster 
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polarization,  about  2  eV,  which  reduces  the  repulsive  contribution  to  the  bonding. 
These  two  contributions  alone  account  for  about  90%  of  the  whole  binding  energy. 
The  other  CSOV  steps,  X  polarization  and  covalency,  only  slightly  improve 
(Table  2).  Similar  data  have  been  obtained  for  Ag4Cl  and  Ag4Br  [5].  The  main 
difference  is  that,  because  of  the  larger  equilibrium  bond  distances,  the  electrostatic 
interaction  and  the  cluster  polarization  are  smaller  than  in  Ag4F.  On  the  other  hand, 
the  contribution  from  the  Cl  and  Br“  polarization  is  higher,  as  expected  on  the  basis 
of  the  larger  dipole  polarizabilities  of  these  anions  (Table  1). 

The  ionic  nature  of  the  silver-halogen  interaction  suggests  that  an  external  electric 
field  F  will  produce  large  changes  on  the  chemisorption  properties.  For  a  cluster  in  a 
uniform  electric  field  F,  the  SCF  first  order  perturbation  theory  energy,  Ep(F),  is 
obtained  as  the  difference  between  the  F  =  0  SCF  energy  and  fix  F,  where  fi  is  the 
field-free  dipole  moment  [12,13].  This  is  a  pure  Stark  effect  and  does  not  include 
chemical  changes  caused  by  the  field.  These  electronic  effects  are  explicitly  considered 
in  the  SCF  variational  energy  in  the  presence  of  the  field,  Escp(F),  obtained  by 
addition  of  the  term  (Sq  —  IQjv;Rm)  x  F  to  the  field-free  hamiltonian;  i  denotes 
elecuons,  N  nuclei,  aind  Q^-  the  nuclear  charges  [12,13]. 

The  fields  considered  are  ±  0.01  a.u.  =  ±  5.7  x  lO^V/cm.  The  sign  of  F  is  such 
that  F  <  0  attracts  electrons  from  the  surface  towards  the  ligand.  Indeed,  we  observe 
that  for  F  =  —  0.01  the  X  ions  are  pulled  away  from  the  surface  giving  a  large 
increase  in  r,  similarly,  F  =  -f-  0.01  pushes  the  ion  toward  the  surface  (Table  3). 
Consequently,  the  vibrational  energies  decrease  for  F  =  -(-  0.01  and  considerably 
increase  for  F  =  —  0.01  (Table  4).  This  occurs  because  the  adsorbed  ligand  is  stretched 
against  the  wall  represented  by  surface  charge  density.  The  changes  caused  by  the 
external  field  are  dominantly  stark  effects,  while  the  chemical  changes  make  only  minor 
contributions  (Tables  3  and  4). 
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In  conclusion,  we  have  shown  that  the  bonding  of  F,  Cl,  and  Br  with  Ag  clusters 
is  almost  entirely  ionic.  We  believe  this  is  true  also  for  surfaces  in  particular 
considering  that  the  SCF  wave  function  underestimates  the  electron  affmities  of  the 
atoms  and  that  a  silver  surface  is  more  polarizable  than  a  small  cluster. 
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Table  1 

Properties  of  Ag4-X  Clusters 


Agj-F 

Ag4-CI 

Ag4-Br 

EA(X) 

eV  SCF 

1.35 

2.08 

2.52 

exp 

3.40 

3.62 

3.36 

«n 

SCF 

1.35 

3.77 

6.19 

HFL* 

1.71 

5.89 

8.18 

r(Ag-X) 

au 

3.310 

4.347 

4.721 

cm“* 

277 

180 

115 

De(Ag4+X) 

eV 

2.72 

2.68 

2.76 

iti(r)  (Ag4-X)  au 

Mo 

-  1.09 

-  1.72 

-2.00 

Ml 

-  1.39 

-  1.60 

-  1.59 

M2 

-  1.09 

-0.04 

-0.01 

M3 

1.01 

0.03 

-0.03 

uir)  fX)  au 

Mo 

-2.90 

-3.31 

-3.43 

Ml 

M, 

-  1.00 

-  1.10 

-  1.11 

0.10 

-0.03 

-0.04 

m; 

-0.01 

0.04 

-0.01 

HFL  =  Hartree-Fock  Limit  [14] 


8 


Table  2 

CSOV  Analysis  for  Ag4-F  and  Ag4-PC  (PC  =  Point  charge)^ 


Step 

p. 

*-int 

A  Hint 

A/i 

0. 

Frozen  orbital 

PC 

F“ 

-r4.82 
+  3.37 

... 

1  1 

1. 

Ag^  polarization  (Ag4  basis) 

PC 

F" 

+  6.92 
+  5.25 

+  2.10 
+  1.88 

-  1.52 

-  1.37 

+  2.97 
+  2.84 

2. 

Ag^  polarization  (full  basis) 

F“ 

+  5.35 

+  0.10 

-1.43 

-0.06 

3. 

X  polarization  (X  basis) 

F~ 

+  5.59 

+  0.24 

-  1.15 

+  0.28 

4. 

X  polarization  (full  basis) 

F” 

+  5.84 

+  0.25 

-  1.05 

+  0.10 

5. 

Covalent  e  sym.  interaction 

F“ 

+  5.85 

+  0.01 

-  1.04 

+  0.01 

6. 

Full  SCF 

F“ 

+  5.88 

+  0.03 

-0.94 

+  0.10 

■"Interaction  energies  in  eV  are  computed  with  respect  to  separate  ions;  in  au 
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Table  3 

Equilibrium  Bond  Distances  in  Bohrs  for  Halogens  Chemisorbed  on  Ag4. 
When  a  Field  F=— 0.01  au  is  Applied  Both  SCF  and  Stark  Values 
Are  Given.  Ar  are  Also  Given 


F 

Ag4-F 

Ag4-Cl 

Ag4-Br 

0.00 

SCF 

3.310/... 

4.347/... 

4.721/... 

+  0.01 

Stark 

3.09S-0.21 

4.055  -0.29 

4.405/ -0.32 

SCF 

3.092  -0.22 

4.05S  -0.29 

4.4 18/+',  30 

-0.01 

Stark 

3.63:’'-0.33 

4.S93  +0.55 

5.375/ +  0.65 

SCF 

3.620' -0.31 

4.904  ^0.56 

5.410/  +  0.69 
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